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Abstract 

The singlc-degcncratc model for the progenitors of type la supernovae (SNc la) is one of the two 
most popular models, in which a carbon-oxygen white dwarf (CO WD) accretes hydrogen-rich material 
from its companion, increases its mass to the Chandrasekhar mass limit, and then explodes as a SN la. 
Incorporating the prescription of Hachisu et al. (1999a) for the accretion efficiency into Eggleton's stellar 
evolution code and assuming that the prescription is valid for all metalhcities, we carried out a detailed 
binary evolution study with different metalhcities. We show the initial and final parameter space for SNe 
la in a (logP — M2) plane. The positions of some famous recurrent novae in the (logP — M2) plane, as well 
as a supersoft X-ray source (SSS), RX J0513. 9-6951 are well explained by our model, and our model can 
also explain the space velocity and mass of Tycho G, which is now suggested to be the companion star of 
Tycho's supernova . Our study indicates that the SSS, V Sge, is a potential progenitor of supernovae like 
SN 2002ic if the delayed dynamical-instability model in Han & Podsiadlowski (2006) is appropriate. 

Key words: stars: binaries: general — stars: supernovae: general — individual(SN 2002ic) — stars: white 
dwarfs 
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1. Introduction 

Type la supernovae (SNe la) play an important role as 
cosmological distance indicators and are successfully ap- 
plied to determine cosmological parameters (e.g. il and 
A; Riess et al. 1998; Perlmutter et al. 1999), which have 
led to the discovery of the accelerating expansion of the 
universe. However, the exact nature of SNe la progenitors 
has not been clarified (see the reviews by Hillebrandt & 
Niemeyer 2000; Leibundgut 2000). It is generally agreed 
that SNe la originate from the thermonuclear runaway 
of a carbon-oxygen white dwarf (CO WD) in a binary 
system. The CO WD accretes material from its compan- 
ion, increases mass to its maximum stable mass, and then 
explodes as a thermonuclear runaway. Almost half of the 
WD mass is converted into radioactive nickel-56 in the ex- 
plosion (Branch 2004), and the amount of nickel-56 deter- 
mines the maximum luminosity of SNe la (Arnett 1982). 

According to the nature of the companions of the mass 
accreting white dwarfs, two competing scenarios have 
been proposed, i.e the double-degenerate channel (DD, 
Iben & Tutukov 1984; Whelan & Iben 1987) and the 
single degenerate channel (SD, Whelan & Iben 1973; 
Nomoto, Thielemann & Yokoi 1984). In the DD chan- 
nel, two CO WDs with a total mass larger than the 
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Chandrasekhar mass limit may coalesce, and then explode 
as a SN la. Although the DD channel is theoretically 
less favored, e.g. double WD mergers are more likely to 
result in accretion-induced collapses rather than SNe la 
(Hillebrandt & Niemeyer 2000), it is premature to exclude 
the channel at present since there exists evidence that the 
channel may contribute to a few SNe la (Howell et al. 
2006; Branch 2006; Quimby, Hofiich & Wheeler 2007). 
The single-degenerate Chandrasekhar model is an alter- 
native of DD and some surveys showed that many SD 
systems are good candidates for SNe la (Whelan & Iben 
1973; Nomoto, Thielemann & Yokoi 1984; Parthasarathy 
et al. 2007). In this model, the maximum stable mass of a 
CO WD is ~ I.378M0 (close to the Chandrasekhar mass, 
Nomoto, Thielemann & Yokoi 1984), and the compan- 
ion is probably a main sequence star or a slightly evolved 
star (WD-I-MS), or a red-giant star (WD-I-RG) (Yungelson 
et al. 1995; Li & van den Heuvel 1997; Hachisu et al. 
1999a; Hachisu et al. 1999b; Nomoto et al. 1999; Langer 
et al. 2000; Han & Podsiadlowski 2004; Meng, Chen & 
Han 2009). The SD model is supported by many obser- 
vations. For example, variable circumstellar absorption 
lines were observed in the spectra of SN la 2006X (Patat 
et al. 2007), which indicates the SD nature of its precur- 
sor. Patat et al. (2007) suggested that the progenitor of 
SN 2006X is a WD + RG system based on the expansion 
velocity of the circumstellar material, while Hachisu et al. 
(2008) argued a WD + MS nature for this SN la. Recently, 
Voss & Nelemans (2008) suggested that SN 2007on is also 
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possible from a WD + MS channel. In this paper, we only 
focus on the WD + MS channel, which is a very impor- 
tant channel for producing SNe la in our Galaxy (Han & 
Podsiadlowski 2004). 

Observationally, some WD + MS systems are possible 
progenitors of SNe la (see the review of Parthasarathy et 
al. 2007). For example, supersoft X-ray sources (SSSs) 
were suggested as good candidates for the progenitors of 
SNe la (Hachisu & Kato 2003b; Hachisu & Kato 2003c). 
Some of the SSSs are WD -I- MS systems and some are 
WD 4- RG systems (Di Stcfano & Kong 2003). Di Stefano 
& Kong (2003) reported that in every galaxy, there are at 
least several hundred SSSs with a luminosity of > 10'^'' erg 
based on Chandra data from four external galaxies (an 
elliptical galaxy, NGC 4967; two face-on spiral galaxies, 
MlOl and M83; and an interacting galaxy, M51). They 
also noticed that the SSSs appear to be associated with the 
spiral arms in the spiral galaxies, which may indicate that 
SSSs are young systems (WD -I- MS systems?). Recurrent 
novae may also be good candidates as the progenitors of 
SNe la (Branch ct al. 1995) and several novae have been 
suggested to be possible progenitors (Hachisu et al. 2000a; 
Hachisu et al. 2000b; Hachisu & Kato 2000; Hachisu & 
Kato 2003a; Hachisu & Kato 2005; Hachisu & Kato 2006a; 
Hachisu & Kato 2006b; Hachisu et al. 2007). Some of the 
novae are WD -t-MS systems and some are WD + RG 
systems. 

A direct way to confirm the progenitor model is to 
search for the companion stars of SNe la in their remnants. 
The discovery of the potential companion of Tycho's su- 
pernova have verified the power of the method and also 
the reliability of the WD + MS model (Tycho G named 
in Ruiz-Lapuente et al. 2004). Recently, Hernandez et al. 
(2009) stressed further the statement in Ruiz-Lapuente ct 
al. (2004) by analysing the chemical abundances of Tycho 
G. 

Many works have concentrated on the WD+MS chan- 
nel. Some authors (Hachisu et al. 1999a; Hachisu et al. 
1999b; Hachisu et al. 2008; Nomoto et al. 1999; Nomoto et 
al. 2003) have studied the WD-I-MS channel by a simple 
analytical method to treat binary interactions. Such an- 
alytic prescriptions may not describe some mass-transfer 
phases, especially those occurring on a thermal time-scale 
(Langer et al. 2000,). Li & van den Heuvel (1997) studied 
this channel from detailed binary evolution calculation, 
while they considered two WD masses (1.0 and 1.2 Mq). 
Langer et al. (2000) investigated the channel for metallic- 
ities Z — 0.001 and 0.02, but they only studied the case 
A evolution (mass transfer during core hydrogen burning 
phase). Han & Podsiadlowski (2004) carried out a de- 
tailed study of this channel including case A and early 
case B (mass transfer occurs at Hertzsprung gap, HG) for 
Z = 0.02. Following the study of Han & Podsiadlowski 
(2004), Meng, Chen & Han (2009) studied the WD + 
MS channel comprehensively and systematically at var- 
ious Z and showed the initial parameter spaces for the 
progenitors of SNe la and the distributions of the initial 
parameters for the progenitors of SNe la. Here, based on 
the study in Meng, Chen & Han (2009), we want to show 
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the final parameter spaces of companions at the moment 
of SNe la explosion and check whether the model used in 
Meng, Chen & Han (2009) can explain the properties of 
some recurrent novae and SSSs, which are suggested to 
be the possible progenitors of SNe la, and the properties 
of Tycho G, which is the potential companion of Tycho's 
supernova. 

In section 2, we simply describe the numerical code for 
binary evolution calculations. The evolutionary results 
are shown in section 3. In section 4, we briefly discuss 
our results, and finally we summarize the main results in 
section 5. 

2. Binary Evolution Calculation 

Our method for treating the binary evolution of WD + 
MS systems is same to that in Meng, Chen & Han (2009). 
In the following, we simply redescribe our method. We 
use the stellar evolution code of Eggleton (1971); Eggleton 
(1972); Eggleton (1973) to calculate the binary evolutions 
of WD-I-MS systems. The code has been updated with 
the latest input physics over the last three decades (Han, 
Podsiadlowski & Eggleton 1994; Pols et al. 1995; Pols 
et al. 1998). Roche lobe overflow (RLOF) is treated 
within the code described by Han et al. (2000). We set 
the ratio of mixing length to local pressure scale height, 
a = I /Hp, to 2.0, and set the convective overshooting 
parameter, Sqv, to 0.12 (Pols et al. 1997; Schroder et 
al. 1997), which roughly corresponds to an overshooting 
length of 0.25i?p. Ten metallicities are adopted here (i.e. 
Z=0.0001, 0.0003, 0.001, 0.004, 0.01, 0.02, 0.03, 0.04, 0.05 
and 0.06). The opacity tables for these metallicties are 
compiled by Chen & Tout (2007) from Iglesias & Rogers 
(1996) and Alexander & Ferguson (1994). For each Z, the 
initial hydrogen mass fraction is obtained from 

X = 0.76-3.0Z, (1) 

since this relation can well reproduce the color-magnitude 
diagram (CMD) of some clusters(Pols et al. 1998). 

Instead of solving stellar structure equations of a WD, 
we adopt the prescription of Hachisu et al. (1999a) on 
WDs accreting hydrogen-rich material from their com- 
panions. In a WD -|- MS system, the companion fills its 
Roche lobe at MS or during HG, and transfers material 
onto the WD. We assume that if the mass-transfer rate, 
IM2I, exceeds a critical value, Mcr, the accreted hydrogen- 
rich material steadily burns on the surface of WD, and is 
converted into helium at the rate of Mcr- The unprocessed 
matter is lost from the system as an optically thick wind 
at a rate of Mwind = IM2I - Mcr (Hachisu et al. 1996). 
Based on the opacity from Iglesias & Rogers (1996), the 
optically thick wind is very sensitive to Fe abundance, 
and it is possible for the wind not to work when Z < 0.002 
(Kobayashi et al. 1998). Thus, there should be an obvious 
low-metallicity threshold for SNe la in comparison with 
SN II. However, this metallicity threshold was not found 
(Prieto et al. 2008). Considering the uncertainties in the 
opacities, we therefore assume rather arbitrarily that the 
optically thick wind is valid for all metallicities. 
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The critical accretion rate is given by 

Afcr = 5.3 X 10~'^ ^^''^~^\ mwu - 0.4), (2) 

where X is hydrogen mass fraction and Mwd is the mass 
of the accreting WD (mass is in Mq and mass-accretion 
rate is in Mq/jt, Hachisu et al. 1999a). The effect of 
metahicities on equation (2) has not been included since 
the effect is very small and can be neglected (Mcng et al. 
2006). 

We adopted the following assumptions when \Al2\ < 
Mcv (1) When IM2I > |Mcr, the hydrogen-sheU burn- 
ing is steady and no mass is lost from the system. (2) 
When iMcr > IM2I > |Mcn a very weak shell flash is trig- 
gered but no mass is lost from the system. (3) When 
IM2I < gAfcr , the hydrogen-shell flash is so strong that 
no material can be accumulated by the accreting CO WD. 
Then, the growth rate of the mass of the helium layer un- 
der the hydrogen-burning shell can be defined as 

A'fHc-?/H|M2|, (3) 

where ?7h is the mass accumulation efficiency for hydrogen 
burning and its values is 

r M„/|Af2|, |M2|>Afer, 
7?H=< 1, A'/cr > |Af2| > |Mer, (4) 

[ 0, |Af2|<|Afer. 

When a certain amount of helium is accumulated, he- 
lium is ignited as a He-flash, and some of the helium is 
blown off from the surface of the CO WD. Then, the mass 
growth rate of the CO WD, Mwd, is 

A/wD = mcMuc = ?7He'^H | Af2 | , (5) 

where t/hc is the mass accumulation efficiency for helium- 
shell flashes, and its value is taken from Kato & Hachisu 
(2004). We assume that if an accreting CO WD increases 
its mass to 1.378 Mq (Nomoto, Thielemann & Yokoi 
1984), it explodes as a SN la. 

3. Results 

3.1. Final State of Binary Evolution 

In figure 1, we show the initial contour for SNe la 
and the final state of binary evolution in (logP — Af2) 
plane at the moment of SNe la explosion for the case of 
M^j3 = 1.1A/q with Z = 0.02. The accreting WD may 
reach a mass of 1.378 Mq during the optically thick wind 
phase (the filled squares. Case Wind) or after optically 
thick wind while in stable (the filled circles. Case Calm) 
or unstable (the filled triangles. Case Nova) hydrogen- 
burning phase. For the first glance, the position of the 
final contour is much lower than that of the initial con- 
tour, which results from the mass transfer from secondary 
to white dwarf and the mass loss from the system. A SSS, 
RX J0513.9-6951 (open star), whose WD mass is 1.30 Mq 
(Hachisu & Kato 2003b), locates in the initial contour. It 
should be a good candidate for the progenitor of SN la, as 
suggested in Hachisu & Kato (2001) and Hachisu & Kato 
(2003b). Three recurrent novae (U Sco, V394 CrA and 
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Fig. 1. Parameter regions producing SNe la in the 
(logP — M2) (orbital period-donor mass) plane for the WD 
-I- MS systems. The initial WD mass is I.IOMq. The 
WD -|- MS system insides the region encircled by sohd hne 
(labeled "initial") will increase its white dwarf mass up to 
1.378 Mq, where we assume a SN la explosion. The fi- 
nal state of the WD -|- MS system in the plane is encir- 
cled by dot-dashed line (labeled "final"). Filled squares indi- 
cate SN la explosions during an optically thick wind phase 
(IM2I > Mcr). Filled circles denote SN la explosions af- 
ter the wind phase, where hydrogen-shell burning is stable 
(Mcr > IM2I > |Mcr). Filled triangles denote SN la explo- 
sions after the wind phase where hydrogen-shell burning is 
mildly unstable {■^Mcr > |Af2| > gMcr). A supersoft X-ray 
sources, RX J0513. 9-6951 (open star) is plotted, whose or- 
bital period is 0.763 days (PakuU et al. 1993). Three recur- 
rent novae are indicated by a filled star for U Sco whose pe- 
riod is 1.2306 days (Schacfer & Ringwald 1995; Hachisu et al. 
2000a; Hachisu et al. 2000b), a solar symbol for CI Aq whose 
orbital period is 0.6184 days (Mennickent & Honcycutt 1995) 
and an earth symbol for V394 CrA whose orbital period is 
0.7577 days (Schaefer 1990). 
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Fig. 2. Similar to Fig. 1, but for different initial WD mass. 

CI Aql) are outside the initial contour for SNe la, while 
still inside the final region. The WD mass of U Sco was 
estimated to be 1.37 Mq (Hachisu et al. 2000a; Hachisu 
et al. 2000b), which is very close to 1.378 Mq. The MS 
mass of U Sco is 1.5 A^© (Hachisu et al. 2000a; Hachisu 
et al. 2000b), which still has enough material to trans- 
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Fig. 3. Similar to Fig. 1, but for different metallicities. 

fer onto the WD and increases the WD mass to 1.378 
Mq. So. U Sco is a very likely candidate for the pro- 
genitor of SN la (sec also Hachisu ct al. 2000a; Hachisu 
et al. 2000b). The WD masses of V394 CrA and CI 
Aql were estimated to be 1.37 Mq and 1.2 Mq, respec- 
tively. Their companion masses are still unclear. The 
best-fit companion masses for the two recurrent novae are 
1.50 Mq, while I.OMq ~ 2.OM0 is still accepted (Hachisu 
& Kato 2000; Hachisu & Kato 2003a; Hachisu, Kato & 
Schaefer 2003). So, based on our binary evolution calcu- 
lation, V394 CrA is a very likely progenitor of SN la and 
CI Aql is a possible progenitor of SN la. A further obser- 
vation to confirm their companion masses is necessary for 
judging their fates. Another very famous recurrent nova is 
V 1974 Cygni (Nova Cygni 1992), whose companion mass 
is 0.21 Mq (Deyoung & Schmidt 1994; Retter, Leibowitz 
& Ofek 1997) and orbital period is 0.0813 days (Paresce 
et al. 1995; Retter, Leibowitz & Ofek 1997). Its position 
is outside the final contour in this paper. Considering its 
low WD mass (I.OMq - I.IMq, Hachisu & Kato 2005), it 
is very unlikely for V 1974 Cygni to become a SN la (See 
Hachisu et al. (2008) for similar discussions about these 
recurrent novae). 

In figure 2, we present the initial contours for SNc la 
and the final state in (log P — M2 ) plane at the moment 
of SNe la explosion for different initial WD masses with 
Z = 0.02. This figure further confirms the likelihood for 
V394 CrA and CI Aq to be the progenitors of SNe la. 

In figure 3, we show the trend of the initial and final 
contour with metallicity. For clarity, we only show the 
case of M^Y) = I.IOMq, since the other cases give similar 
results (see Fig. 4 in Meng, Chen & Han 2009). We see 
from the figure that there is a clear trend for the initial 
and final contour to move to higher mass with metallicity. 
This is due to the correlation between stellar evolution and 
metallicity. Generally, stars with high metallicity evolve 
in a way similar to those with low metallicity but less mass 
(Umcda et al. 1999; Chen & Tout 2007; Meng, Chen & 
Han 2008). Thus, for binaries of CO WDs with partic- 
ular orbital periods, the companion mass increases with 
metallicity (see Meng, Chen & Han 2009 in details). 
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Fig. 4. Similar to figure f, but only for Case Wind. The 
initial WD mass is f.fO Mq. The big field squares show 
the initial WD -|- MS systems, while the small ones show 
their final positions at the moment of SNe la explosion. 
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Fig. 5. Similar to figure 4, but for different initial WD 
mass. Different colors represent different initial WD masses. 
The big field squares show the initial WD -|- MS systems, 
while the small ones show their final positions at the mo- 
ment of SNe la explosion. In the figure, we do not plot 
the cases with Mi < O.QMq for the clarity of appearance 
of the figure since no Case Wind occurs for the cases. 

We present the initial and final state of binary evolution 
of SNe la for Case Wind in figures 4 and 5. Figure 4 is only 
for the case of M^^} = I-IMq and figure 5 is for different 
initial WD masses. Here, we only show the case oi Z = 0.02 
as a typifier, since the other cases give similar results (see 
figure 3). In these two figure, we see that the range of 
initial MS masses with Z = 0.02 for the Case Wind is be- 
tween 2.2 Mq and 3.5 Mq. However, the range strongly 
depends on metallicity, that is: the progenitor systems 
have more massive companions for a higher Z (see figure 
3 in this paper and Fig. 4 in Meng, Chen & Han 2009). 
Hachisu et al. (2008) found that when M^j^^ < 0.9Mq, no 
Case Wind exists for Z = 0.02. We get a similar results 
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for Z = 0.02. However, this result also rests on metallic- 
ity. For example, the mass limit for no wind is 0.8 A/© for 
Z = 0.06, while it is larger than 1.0 Mq for Z = 0.0001. 
In addition, Hachisu et al. (2008) noticed that only when 
M\ > S.OMq, WDs explode at the optically thick wind 
phase as SNe la. Their result is much different from that 
in this paper. The difference is derived from the differ- 
ent treatment on mass transfer between WD and MS and 
from different mass-loss mechanism from systems. Firstly, 
they used a analytic method for treating binary interac- 
tion. This method can not describe certain mass-transfer 
phase, in particular those occurring on a thermal time- 
scale and may overestimate mass-transfer phase on the 
thermal time-scale (Langer et al. 2000; see also the Figs. 
1 and 4 in Han & Podsiadlowski 2004). Secondly, they as- 
sume a mass-stripping effect, that is: optically thick winds 
from WD collide with secondary and strip off its surface 
layer. This effect may attenuate the mass-transfer rate be- 
tween WD and its companion, but the total mass-loss rate 
from secondary increases (Hachisu & Kato 2003a; Hachisu 
& Kato 2003b; Hachisu et al. 2008). For the two reasons 
above, more materials thus are lost from binary systems 
as the optically thick wind and stripped-off material be- 
fore supernova explosion, and then the systems are more 
likely to explode after the optically thick wind. 

The material lost as the optically think wind forms cir- 
cumstellar material (CSM) which may result in a color 
excess of SNe la (Meng et al. 2009). Since the CSM is 
very close to SN la for Case Wind, the ejecta of explosion 
may interact with the CSM, and the SN la might then be 
observed as supernovae like SN 2006X which show a vari- 
able Na I D line (Patat et al. 2007; Blondin et al. 2009). 
We will discuss this in the futiu'c paper. 

In the final contours, there is no special region for Case 
Wind (see figure 5), that is: WDs can explode as SNe 
la at any position in the final permitted region in the 
(logP — M2) plane, which is much different from the Case 
Calm and the Case Nova (see the following subsection). 

3.3. Case Calm 

If the mass-transfer rate from secondary onto WD is 
below the critical accretion rate (Equation 2) while above 
the lowest rate of steady hydrogen burning before super- 
nova explosion, i.e. M^r > IM2I > \Mcr, the WD under- 
goes steady H burning at the moment of SN la explosion. 
This is the reason why we call this "Case Calm" (Same to 
that in Hachisu et al. 2008). During the steady hydrogen 
burning phase, the WD may be observed as a SSS. The 
material lost as the optically thick wind forms CSM, but it 
has been dispersed too washy to be detected immediately 
after the SN la explosion. In figures 6 and 7, we show 
the initial and final state of binary evolution of SNe la for 
Case Calm in the (logP — M2) plane. Figure 6 is only for 
the case of = 1.1 Mq and figure 7 is for different ini- 

tial WD masses. Here, we only show the case oi Z = 0.02 
as a typifier, since the other cases give similar results (see 
also figure 3). In figure 6, it is clear that there is gap 
for both initial system and final state in the (logP — M2) 
plane, and the gap divides the Case Calm into two groups. 
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Fig. 6. Similar to figure 4, but only for Case 
Calm. The big field circles show the initial WD 
-I- MS systems, while the small ones show their fi- 
nal positions at the moment of SNe la explosion. 
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Fig. 7. Similar to Fig. 6, but for different initial 
WD masses. Different colors represent different ini- 
tial WD masses. The big field circles show the ini- 
tial WD -|- MS systems, while the small ones show 
their final positions at the moment of SNe la explosion. 

The gap in the initial contour is the region for Case wind 
and the area of the gap decrease with decreasing. 
When M^j-, < O.9M0, the gap disappear (see figure 7). 

The gap is also shown in the initial (logP — M2) plane 
of Han & Podsiadlowski (2004) and Meng, Chen & Han 
(2009), but they do not discussed it. In this paper, we 
explain its origin as foUowings. For the high-mass group 
in the initial contour of figure 6, their mass ratio is large. 
After the onset of Roche lobe overflow (RLOV), mass- 
transfer rate is so high that the system almost undergoes 
dynamically unstable. At this stage, a large amount of 
hydrogen-rich materials lose as the optically thick wind 
and hence the secondary mass decreases sharply. The 
mass-transfer rate drops after the mass ratio has been 
reversed and then the optically thick wind stops. When 
CO WD mass increases to 1.378 Mq, the secondary mass 
decreases to about 0.9 A/q, which corresponds to the low- 
mass group in the final state contour of figure 6. 
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For the low-mass group in the initial contour of figure 6, 
their mass ratio is not very large and only a small amount 
of hydrogen-rich materials lose as the optically thick wind, 
which means that most of transferred materials are accu- 
mulated on the WD. Then, the secondary mass decreases 
slightly, which forms the high-mass group in the final state 
contour of figure 6. 

In figure 7, there is a blank region in the final state re- 
gion for Case Calm (around Ml^ = 1.4A/Q-logP = 0.6), 
which means that SSSs may not be observed at this re- 
gion before SNe la explosion. This phenomena is much 
different from that of Case Wind. The blank region is 
mainly derived from the low-mass group in the initial con- 
tour. For a system belonging to the group, mass transfer 
is almost conservative, and then the period is always de- 
creasing since mass ratio is not reversed, which leads to 
the region absenting long-period systems. For a system 
in the high-mass, although its period decrease during op- 
tically thick wind, the period increases after mass-ratio 
reversion since no material loses from system and mass 
transfer becomes conservative (see the panel (h) of Fig.l 
in Han & Podsiadlowski 2004). Thus, its final period is a 
bit complicated. 

34- Case Nova 



M\,, 
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Fig. 8. Similar to figure 4, but only for Case Nova. 
The big field triangles show the initial WD -|- 
MS systems, while the small ones show their fi- 
nal positions at the moment of SNe la explosion. 



If mass-transfer rate from secondary onto WD is be- 
low the lowest rate of steady hydrogen burning, while still 
higher than a special value, i.e. ■^Mcr > IM2I > |-^cr, 
hydrogen shell burning is unstable to flash, and this can 
recur many times in a short period as a recurrent nova. 
We call this "Case Nova" (similar to Hachisu et al. 2008). 
During the last recurrent nova phase, the amount of the 
hydrogen-rich material lost from system is very small. On 
the other hand, the material lost as the optically thick 
wind is already much far from the center of SN la explo- 
sion. It will take 100-1000 yr for the SN la ejecta to catch 
the material lost as the optically thick wind (Hachisu et 
al. 2008), and then no signal about CSM can be detected 




-1 -0.5 0.5 1 

logP(dB7a) 

Fig. 9. Similar to figure 8, but for different initial 
WD masses. Different colors represent different ini- 
tial WD masses. The big field circles show the ini- 
tial WD + MS systems, while the small ones show 
their final positions at the moment of SNe la explosion. 



immediately after supernova explosion. In figures 8 and 
9, we show the initial and final state of binary evolution 
of SNe la for Case Nova in the (logP — M2) plane. Figure 

8 is only for the case of M^^) = l-l-^o and figure 9 is 
for different initial WD masses. Here, we only show the 
case oi Z = 0.02 as a typifier, since the other cases give 
similar results (see also figure 3). The initial region of 
Case Nova is also divided into two group, which is similar 
to that of Case Calm. The systems located in the gap 
are those for Case Wind or Case Calm. The difference 
between Case Calm and Case Nova is that the gap dis- 
appear when = O.8M0, not 0.9Mo (see figure 9). 
For a similar reason to that of Case Calm, the high-mass 
group in the initial contour forms the low-mass group in 
the final contour, while the low-mass group in the initial 
contour forms the high-mass group in the final contour. 

Similar to figure 7, there is also a blank region in figure 

9 (around Aff^ = 1.3AfQ-logP = 0.5). The origin of the 
blank region is similar to that of Case Calm. 

The recurrent nova U Sco, which is an excellent candi- 
date of SN la progenitor (Hachisu et al. 2000a; Hachisu et 
al. 2000b), is located in the middle of the final region, and 
its position in the (logP — M2) plane is just the permitted 
region for recurrent nova (see figure 9). Our model can 
explain the position of U Sco excellently. 

3.5. Companion state after SNe la explosion. 

A good way of discriminating between the many SN la 
progenitor scenarios is to search for the companion of a 
SN la in its remnant. Unless the companion is another 
WD (DD channel; in which it has been destroyed by the 
mass-transfer process itself before explosion), it survives 
and shows some special properties in its spectra, which 
is originated from the contamination of supernova ejecta 
(Marietta et al. 2000; Ruiz-Lapuente et al. 2004; Branch 
2004). Tycho's supernova, which is one of only two SNe 
la observed in our Galaxy, provides an opportunity to 
address observationally the identification of the surviving 
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Fig. 10. The final state of companion star at the mo- 
ment of SNc la explosion for different initial WD mass 
with Z = 0.02. Cross represents Tycho G, which 
is a potential candidate for the companion of Tycho's 
supernova (Ruiz-Lapuente et al. 2004; Branch 2004). 
The length of the cross represents observational error. 
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companion. Ruiz-Lapuente et al. (2004) searched the re- 
gion of the remnant of Tycho's supernova and suggested 
that Tycho G, a sun-hke star, is the companion of Tycho's 
supernova. Although, Ihara et al. (2007) argued that the 
spectrum of Tycho G do not show any special properties, 
which seems to exclude the possibility of Tycho G to be 
the companion of Tycho' supernova, the analysis of the 
chemical abundances of the Tycho G upholds the com- 
panion nature of the Tycho G (Hernandez et al. 2009). 

However, the knowledge about the companions of SNe 
la after explosions is still unclear. Generally, the super- 
nova ejecta in the single degenerate model collides into the 
envelope of its companion and strips some hydrogen-rich 
material from the surface of the companion. After the col- 
lision, the companion gains a kick velocity, which is much 
smaller than orbital velocity, and leaves explosion center 
at a velocity similar to its orbital velocity (Marietta et al. 
2000; Meng et al. 2007). 

For a companion star in the remnant of a SN la, its 



mass and space velocity can be detected directly. Since 
the stripped-off hydrogen-rich material from companion 
surface and the change of the space velocity of companion 
resulting from the collision of explosion ejecta are both 
small (Marietta et al. 2000; Mattila et al. 2005; Leonard 
2007; Meng et al. 2007), we can approximately use the 
state of a companion at the moment of supernova ex- 
plosion to present its final state after supernova explo- 
sion. In figure 10, the final states of companion stars are 
presented in Voth ~ Mf (orbital velocity-final companion 
mass) plane for different initial WD mass with Z = 0.02. 
From the figure, we can see that the range of space veloc- 
ity is from 70 km/s to 210 km/s and the mass is between 
0.6 Mq and 2.5 Mq. Tycho G just locates in this range. 
Our work can explain the observation. In the future pa- 
per, we will give detailed results about Tycho G by the 
approach of binary population synthesis. 

The effect of metallicity on the final state of the com- 
panions are shown in figure 11. In the figure, we only 
show the case of M^j-, = I.IOMq, since the other cases 
give similar results, we can see from the figure that all 
the metallicity may explain the properties of Tycho G. 
But for a detailed binary population synthesis study, only 
the case of Z = 0.02 can well explain the position of Tycho 
G in the T4rb — Mf ^ plane. Please notice our following 
paper. 

4. Discussion 

Initial and Final Contour 

In this paper, we show the initial and final parameter 
spaces for SNe la in (logP — M2) plane. The final masses 
of the companions are between 0.6 Mq and 2.5 Mq for 
Z = 0.02. However, the range of the companion mass in 
Hachisu et al. (2008) is from 1.2 Mq to about 3 Mq. 
This difference is directly derived from the different ini- 
tial contours, which is mainly determined by the different 
treatment of thermal timescale mass transfer and different 
mass-loss mechanism. They use a simply analytic method 
to estimate the mass-transfer rate, which leads to a differ- 
ent initial contour for SNe la, especially for low-mass WD 
(Han & Podsiadlowski 2004). In addition, they assumed 
a mass-stripping effect, and then, the initial companion 
mass may be as large as 8 A/q, which depends on the 
efficiency of the mass-stripping effect. Since their final 
mass of companion is always larger than 1.2 A/q, it is dif- 
ficult for their model to explain the properties of Tycho 
G; which remains consistent with the surviving compan- 
ion of Tycho's supernova (1 Mq, Ruiz-Lapuente et al. 
2004; Hernandez et al. 2009). Since our results are based 
on detailed binary evolution calculation with latest input 
physics, our results are more realistic and the final con- 
tours in this paper are more likely to approach the real 
ones than those in Hachisu et al. (2008) (see also in Han 
& Podsiadlowski 2004) . Because of the difference range of 
final companion as mentioned above, Tycho G can be a 
likely progenitor of Tycho's supernova in our model, but 
can not be in the model of Hachisu et al. (2008). 
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4-2. Supernova like SN 2002ic and V Sge 
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Fig. 12. Initial region for SNc la with different metallici- 
ties. Here, only the cases of Af^^ = 1.2OM0 are shown. 
The black solid line shows the initial region of the progen- 
itor for supernovae like SN 2002ic. Two supersoft X-ray 
sources, RX J0513. 9-6951 (open star) and V Sge (filled star), 
are plotted, whose masses are 2.7 Mq and 3.5 Mq (Hachisu 
& Kato 2003b; Hachisu &; Kato 2003c), and whose orbital 
period are 0.76 days and 0.51 days, respectively (PakuU 
et al. 1993; Herbig et al. 1965; Patterson ct al. 1998). 

Until SN 2002ic was discovered (Hamuy et al. 2003), 
it was long believed that there are no hydrogen lines in 
the spectra of SNe la. The strong hydrogen lines in the 
spectra of SN 2002ic were explained by the interaction be- 
tween the SN ejecta and the circumstellar material (CSM) 
(Hamuy et al. 2003). Recently, two co-twins of SN 2002ic 
(SN 2005gj and SN 2006gy) were also found (Aldering 
et al. 2006; Ofck et al. 2007). To explain these rare 
objects, many models were suggested and here we just 
list some of them: Hachisu et al. (1999b); Hamuy et al. 
(2003); Livio & Riess (2003); Chugai & Yungelson (2004); 
Han & Podsiadlowski (2006); Hachisu ct al. (2008); Lii 
et al. (2009). Among all the models listed above, the 
delayed dynamical instability model suggested by Han & 
Podsiadlowski (2006) is more interesting and some pre- 
dictions from the model seem to be consistent with obser- 
vations, especially in a sense of the birth rate and delay 
time of these rare objects (Aldering ct al. 2006; Prieto 
et al. 2007). In the scenario of Han & Podsiadlowski 
(2006), SN 2002ic might be from the WD -f- MS channel, 
where the CO WD accretes material from its relatively 
massive companion ('^ 3.OM0), and increases its mass to 
^ 1.3OA'/0 before experiencing a delayed dynamical insta- 
bility. Adopting their model, Meng, Chen & Han (2009) 
assumed that the CO WD can increase its mass to 1.378 
Mq and explode as a SN 2002ic-like case if the mass of 
the CO WD exceeds 1.30 Mq before the delayed dynami- 
cal instability, and then predicted that supernovae like SN 
2002ic may not be found in extremely low-metallicity envi- 
ronments. Under the same assumption to Meng, Chen & 
Han (2009) , we show the permitted region for the delayed 
dynamical instability model with various metallicities in 
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figure 12. For the simplicity of the assumption, the real 
region for the delayed dynamical instability model may 
be larger than that shown in figure 12 (see Fig. 5 in 
Han & Podsiadlowski 2006 and Fig. 2. in Meng, Chen 
& Han 2009). It is clear that the permitted region for 
supernova like SN 2002ic locates above the contours for 
normal supernova. There is a overlap between supernova 
like SN 2002ic and normal SN la, since a progenitor of 
supernova like SN 2002ic with low metallicity may have 
a lower mass than that of normal supernova. Whatever, 
more evidence is necessary to confirm the scenario in Han 
& Podsiadlowski (2006), especially to find a progenitor 
system at present. 

V Sge is a well observed quasi-periodic transient SSS 
in our Galaxy. It switches optical state between a high 
iy ~ 11 mag) and low {y ^ 12 mag) state, and during 
the low state, very soft and very weak X-ray can be de- 
tected (Greiner & van Tecschng 1998). The mass- loss 
rate for V Sge can be as large as ~ lO~^M0yr~^ indi- 
cated by radio observation. (Lockley et al. 1997; Lockley 
et al. 1999). Hachisu & Kato (2003c) use optically thick 
wind to explain the observational property of V Sge excel- 
lently, and in their model, the mass-loss rate may reaches 
as high as ^ lO~^M0yr~^, which is consistent with radio 
observations. Hachisu & Kato (2003c) also estimated the 
WD mass as Mwd ~ 1.25M0 and its companion mass as 
M2 ^ 3.5M0, and then they suggested that V Sge will ex- 
plode as a SN la after about 10~^yr. However, they must 
assume a mass-stripping effect, otherwise the position of 
V Sge in (logP — M2) plane will beyond the permitted re- 
gion for stable mass transfer (Hachisu et al. 2008). We get 
a similar result to that in Hachisu et al. (2008) that V Sge 
locates in a forbidden zone for stable mass transfer if there 
is no special assumption (see in figure 12). However, it lo- 
cates at the boundary of the contour for supernovae like 
SN 2002ic as indicated in Han & Podsiadlowski (2006). 
Considering the smaller region here than a real one and 
a moderately larger WD mass of V Sge than that in this 
paper, V Sge is a potential progenitor system that will 
undergo the delayed dynamical instability as suggested 
in Han & Podsiadlowski (2006). V Sge thus would not 
explode as a normal SN la, but one like SN 2002ic, if 
the model in Han & Podsiadlowski (2006) is appropriate. 
Then, we might regard "V Sge-type star" as the progen- 
itor of supernovae like SN 2002ic if it belongs to WD + 
MS system. 

The lifetime of V Sge-type star is typically ^ lO^yr, 
which is mainly due to the time-averaged wind mass- 
loss rate of ~ IQ-^Mqyt-'^ (Han & Podsiadlowski 2006). 
Since no more than 1 in 100 SNe la belong to the sub- 
group of 2002ic-like supernova (Han & Podsiadlowski 
2006; Aldering et al. 2006; Prieto et al. 2007), there should 
be several V Sge-type star belonging to WD -I- MS sys- 
tem in our Galaxy if we take the Galactic birth rate of 
SNe la as 3-4x 10~^yr~'^ (van den Bergh & Tammann 
1991; CappcUaro & Turatto 1997). Steiner & Diaz (1998) 
listed four V Sge-type stars in the Galaxy and discussed 
their similar spectroscopic and photometric properties. At 
present, their companion masses have not been identified. 
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while their orbital periods is in the range of 0.2-0.5 days, 
which falls into the orbital periods predicted by our model 
(see the final regions in figure 3). Please notice that not 
all SSSs belong to WD + MS systems and the nature of 
the compact star still remains an open question (Steiner 
& Diaz 1998). White dwarfs, neutron stars, black holes, 
or even carbon main-sequence stars are all possible candi- 
dates (Steiner & Diaz 1998). Another alternative for the 
SSSs is WD -I- He star system, which may also contribute 
to a few part of SNe la (Wang et al. 2009). Thus, consid- 
ering the uncertainty of the nature of the four V Sge-type 
star, our estimation on the number of V Sge-type stars is 
consistent with observations. 

However, although the life, the mass-lose rate and the 
number of V Sge-type star seem to match with the pre- 
diction from the delayed dynamical instability model, the 
case of = 1.2Mq is the only one in our model grids 

to account for the position of V Sge in the (logP — M2) 
plane, while there exist arguments on whether the mass of 
a CO WD may be as large as 1.2Mo and whether the CO 
WD of I.2M0 may explode as a SN la (Umeda et al. 1999; 
Meng, Chen & Han 2008). So, obtaining a conclusion on 
whether there is a relation between V Sge-type star and 
the delayed dynamical instability model is premature, and 
V Sge is only a potential candidate for the progenitor of 
supernovae like SN 2002ic. 

5. Summary and Conclusions 

Adopting the prescription in Hachisu et al. (1999a) for 
the mass accretion of CO WDs, we have carried out de- 
tailed binary evolution calculations for the progenitors of 
SNe la with different metallicity in the single-degenerate 
channel (the WD + MS channel) and obtained the initial 
and final parameters in the (logP — A/2) plane. Our model 
may explain the positions of some famous recurrent novae 
in the (logP - M2) plane, as well as a SSS, RX J0513.9- 
6951. Our model can also explain the space velocity and 
mass of Tycho G, which is now suggested to be the poten- 
tial companion star of Tycho's supernova (Ruiz-Lapuente 
et al. 2004; Hernandez et al. 2009). Based on the de- 
layed dynamically instable model in Han & Podsiadlowski 
(2006), we might regard V Sge as a potential progenitor 
of supernova like SN 2002ic. 
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Abstract 

The singlc-degcncratc model for the progenitors of type la supernovae (SNc la) is one of the two 
most popular models, in which a carbon-oxygen white dwarf (CO WD) accretes hydrogen-rich material 
from its companion, increases its mass to the Chandrasekhar mass limit, and then explodes as a SN la. 
Incorporating the prescription of Hachisu et al. (1999a) for the accretion efficiency into Eggleton's stellar 
evolution code and assuming that the prescription is valid for all metalhcities, we carried out a detailed 
binary evolution study with different metalhcities. We show the initial and final parameter space for SNe 
la in a (logP — M2) plane. The positions of some famous recurrent novae in the (logP — M2) plane, as well 
as a supersoft X-ray source (SSS), RX J0513. 9-6951 are well explained by our model, and our model can 
also explain the space velocity and mass of Tycho G, which is now suggested to be the companion star of 
Tycho's supernova . Our study indicates that the SSS, V Sge, is a potential progenitor of supernovae like 
SN 2002ic if the delayed dynamical-instability model in Han & Podsiadlowski (2006) is appropriate. 

Key words: stars: binaries: general — stars: supernovae: general — individual(SN 2002ic) — stars: white 
dwarfs 



> 

o 

00 

06 
O 

o\ 
o 



1. Introduction 

Type la supernovae (SNe la) play an important role as 
cosmological distance indicators and are successfully ap- 
plied to determine cosmological parameters (e.g. il and 
A; Riess et al. 1998; Perlmutter et al. 1999), which have 
led to the discovery of the accelerating expansion of the 
universe. However, the exact nature of SNe la progenitors 
has not been clarified (see the reviews by Hillebrandt & 
Niemeyer 2000; Leibundgut 2000). It is generally agreed 
that SNe la originate from the thermonuclear runaway 
of a carbon-oxygen white dwarf (CO WD) in a binary 
system. The CO WD accretes material from its compan- 
ion, increases mass to its maximum stable mass, and then 
explodes as a thermonuclear runaway. Almost half of the 
WD mass is converted into radioactive nickel-56 in the ex- 
plosion (Branch 2004), and the amount of nickel-56 deter- 
mines the maximum luminosity of SNe la (Arnett 1982). 

According to the nature of the companions of the mass 
accreting white dwarfs, two competing scenarios have 
been proposed, i.e the double-degenerate channel (DD, 
Iben & Tutukov 1984; Whelan & Iben 1987) and the 
single degenerate channel (SD, Whelan & Iben 1973; 
Nomoto, Thielemann & Yokoi 1984). In the DD chan- 
nel, two CO WDs with a total mass larger than the 
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Chandrasekhar mass limit may coalesce, and then explode 
as a SN la. Although the DD channel is theoretically 
less favored, e.g. double WD mergers are more likely to 
result in accretion-induced collapses rather than SNe la 
(Hillebrandt & Niemeyer 2000), it is premature to exclude 
the channel at present since there exists evidence that the 
channel may contribute to a few SNe la (Howell et al. 
2006; Branch 2006; Quimby, Hofiich & Wheeler 2007). 
The single-degenerate Chandrasekhar model is an alter- 
native of DD and some surveys showed that many SD 
systems are good candidates for SNe la (Whelan & Iben 
1973; Nomoto, Thielemann & Yokoi 1984; Parthasarathy 
et al. 2007). In this model, the maximum stable mass of a 
CO WD is ~ I.378M0 (close to the Chandrasekhar mass, 
Nomoto, Thielemann & Yokoi 1984), and the compan- 
ion is probably a main sequence star or a slightly evolved 
star (WD-I-MS), or a red-giant star (WD-I-RG) (Yungelson 
et al. 1995; Li & van den Heuvel 1997; Hachisu et al. 
1999a; Hachisu et al. 1999b; Nomoto et al. 1999; Langer 
et al. 2000; Han & Podsiadlowski 2004; Meng, Chen & 
Han 2009). The SD model is supported by many obser- 
vations. For example, variable circumstellar absorption 
lines were observed in the spectra of SN la 2006X (Patat 
et al. 2007), which indicates the SD nature of its precur- 
sor. Patat et al. (2007) suggested that the progenitor of 
SN 2006X is a WD + RG system based on the expansion 
velocity of the circumstellar material, while Hachisu et al. 
(2008) argued a WD + MS nature for this SN la. Recently, 
Voss & Nelemans (2008) suggested that SN 2007on is also 



2 

possible from a WD + MS channel. In this paper, we only 
focus on the WD + MS channel, which is a very impor- 
tant channel for producing SNe la in our Galaxy (Han & 
Podsiadlowski 2004). 

Observationally, some WD + MS systems are possible 
progenitors of SNe la (see the review of Parthasarathy et 
al. 2007). For example, supersoft X-ray sources (SSSs) 
were suggested as good candidates for the progenitors of 
SNe la (Hachisu & Kato 2003b; Hachisu & Kato 2003c). 
Some of the SSSs are WD -I- MS systems and some are 
WD 4- RG systems (Di Stcfano & Kong 2003). Di Stefano 
& Kong (2003) reported that in every galaxy, there are at 
least several hundred SSSs with a luminosity of > 10'^'' erg 
based on Chandra data from four external galaxies (an 
elliptical galaxy, NGC 4967; two face-on spiral galaxies, 
MlOl and M83; and an interacting galaxy, M51). They 
also noticed that the SSSs appear to be associated with the 
spiral arms in the spiral galaxies, which may indicate that 
SSSs are young systems (WD -I- MS systems?). Recurrent 
novae may also be good candidates as the progenitors of 
SNe la (Branch ct al. 1995) and several novae have been 
suggested to be possible progenitors (Hachisu et al. 2000a; 
Hachisu et al. 2000b; Hachisu & Kato 2000; Hachisu & 
Kato 2003a; Hachisu & Kato 2005; Hachisu & Kato 2006a; 
Hachisu & Kato 2006b; Hachisu et al. 2007). Some of the 
novae are WD -t-MS systems and some are WD + RG 
systems. 

A direct way to confirm the progenitor model is to 
search for the companion stars of SNe la in their remnants. 
The discovery of the potential companion of Tycho's su- 
pernova have verified the power of the method and also 
the reliability of the WD + MS model (Tycho G named 
in Ruiz-Lapuente et al. 2004). Recently, Hernandez et al. 
(2009) stressed further the statement in Ruiz-Lapuente ct 
al. (2004) by analysing the chemical abundances of Tycho 
G. 

Many works have concentrated on the WD+MS chan- 
nel. Some authors (Hachisu et al. 1999a; Hachisu et al. 
1999b; Hachisu et al. 2008; Nomoto et al. 1999; Nomoto et 
al. 2003) have studied the WD-I-MS channel by a simple 
analytical method to treat binary interactions. Such an- 
alytic prescriptions may not describe some mass-transfer 
phases, especially those occurring on a thermal time-scale 
(Langer et al. 2000,). Li & van den Heuvel (1997) studied 
this channel from detailed binary evolution calculation, 
while they considered two WD masses (1.0 and 1.2 Mq). 
Langer et al. (2000) investigated the channel for metallic- 
ities Z — 0.001 and 0.02, but they only studied the case 
A evolution (mass transfer during core hydrogen burning 
phase). Han & Podsiadlowski (2004) carried out a de- 
tailed study of this channel including case A and early 
case B (mass transfer occurs at Hertzsprung gap, HG) for 
Z = 0.02. Following the study of Han & Podsiadlowski 
(2004), Meng, Chen & Han (2009) studied the WD + 
MS channel comprehensively and systematically at var- 
ious Z and showed the initial parameter spaces for the 
progenitors of SNe la and the distributions of the initial 
parameters for the progenitors of SNe la. Here, based on 
the study in Meng, Chen & Han (2009), we want to show 
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the final parameter spaces of companions at the moment 
of SNe la explosion and check whether the model used in 
Meng, Chen & Han (2009) can explain the properties of 
some recurrent novae and SSSs, which are suggested to 
be the possible progenitors of SNe la, and the properties 
of Tycho G, which is the potential companion of Tycho's 
supernova. 

In section 2, we simply describe the numerical code for 
binary evolution calculations. The evolutionary results 
are shown in section 3. In section 4, we briefly discuss 
our results, and finally we summarize the main results in 
section 5. 

2. Binary Evolution Calculation 

Our method for treating the binary evolution of WD + 
MS systems is same to that in Meng, Chen & Han (2009). 
In the following, we simply redescribe our method. We 
use the stellar evolution code of Eggleton (1971); Eggleton 
(1972); Eggleton (1973) to calculate the binary evolutions 
of WD-I-MS systems. The code has been updated with 
the latest input physics over the last three decades (Han, 
Podsiadlowski & Eggleton 1994; Pols et al. 1995; Pols 
et al. 1998). Roche lobe overflow (RLOF) is treated 
within the code described by Han et al. (2000). We set 
the ratio of mixing length to local pressure scale height, 
a = I /Hp, to 2.0, and set the convective overshooting 
parameter, Sqv, to 0.12 (Pols et al. 1997; Schroder et 
al. 1997), which roughly corresponds to an overshooting 
length of 0.25i?p. Ten metallicities are adopted here (i.e. 
Z=0.0001, 0.0003, 0.001, 0.004, 0.01, 0.02, 0.03, 0.04, 0.05 
and 0.06). The opacity tables for these metallicties are 
compiled by Chen & Tout (2007) from Iglesias & Rogers 
(1996) and Alexander & Ferguson (1994). For each Z, the 
initial hydrogen mass fraction is obtained from 

X = 0.76-3.0Z, (1) 

since this relation can well reproduce the color-magnitude 
diagram (CMD) of some clusters(Pols et al. 1998). 

Instead of solving stellar structure equations of a WD, 
we adopt the prescription of Hachisu et al. (1999a) on 
WDs accreting hydrogen-rich material from their com- 
panions. In a WD -|- MS system, the companion fills its 
Roche lobe at MS or during HG, and transfers material 
onto the WD. We assume that if the mass-transfer rate, 
IM2I, exceeds a critical value, Mcr, the accreted hydrogen- 
rich material steadily burns on the surface of WD, and is 
converted into helium at the rate of Mcr- The unprocessed 
matter is lost from the system as an optically thick wind 
at a rate of Mwind = IM2I - Mcr (Hachisu et al. 1996). 
Based on the opacity from Iglesias & Rogers (1996), the 
optically thick wind is very sensitive to Fe abundance, 
and it is possible for the wind not to work when Z < 0.002 
(Kobayashi et al. 1998). Thus, there should be an obvious 
low-metallicity threshold for SNe la in comparison with 
SN II. However, this metallicity threshold was not found 
(Prieto et al. 2008). Considering the uncertainties in the 
opacities, we therefore assume rather arbitrarily that the 
optically thick wind is valid for all metallicities. 
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The critical accretion rate is given by 

Afcr = 5.3 X 10~'^ ^^''^~^\ mwu - 0.4), (2) 

where X is hydrogen mass fraction and Mwd is the mass 
of the accreting WD (mass is in Mq and mass-accretion 
rate is in Mq/jt, Hachisu et al. 1999a). The effect of 
metahicities on equation (2) has not been included since 
the effect is very small and can be neglected (Mcng et al. 
2006). 

We adopted the following assumptions when \Al2\ < 
Mcv (1) When IM2I > |Mcr, the hydrogen-sheU burn- 
ing is steady and no mass is lost from the system. (2) 
When iMcr > IM2I > |Mcn a very weak shell flash is trig- 
gered but no mass is lost from the system. (3) When 
IM2I < gAfcr , the hydrogen-shell flash is so strong that 
no material can be accumulated by the accreting CO WD. 
Then, the growth rate of the mass of the helium layer un- 
der the hydrogen-burning shell can be defined as 

A'fHc-?/H|M2|, (3) 

where ?7h is the mass accumulation efficiency for hydrogen 
burning and its values is 

r M„/|Af2|, |M2|>Afer, 
7?H=< 1, A'/cr > |Af2| > |Mer, (4) 

[ 0, |Af2|<|Afer. 

When a certain amount of helium is accumulated, he- 
lium is ignited as a He-flash, and some of the helium is 
blown off from the surface of the CO WD. Then, the mass 
growth rate of the CO WD, Mwd, is 

A/wD = mcMuc = ?7He'^H | Af2 | , (5) 

where t/hc is the mass accumulation efficiency for helium- 
shell flashes, and its value is taken from Kato & Hachisu 
(2004). We assume that if an accreting CO WD increases 
its mass to 1.378 Mq (Nomoto, Thielemann & Yokoi 
1984), it explodes as a SN la. 

3. Results 

3.1. Final State of Binary Evolution 

In figure 1, we show the initial contour for SNe la 
and the final state of binary evolution in (logP — Af2) 
plane at the moment of SNe la explosion for the case of 
M^j3 = 1.1A/q with Z = 0.02. The accreting WD may 
reach a mass of 1.378 Mq during the optically thick wind 
phase (the filled squares. Case Wind) or after optically 
thick wind while in stable (the filled circles. Case Calm) 
or unstable (the filled triangles. Case Nova) hydrogen- 
burning phase. For the first glance, the position of the 
final contour is much lower than that of the initial con- 
tour, which results from the mass transfer from secondary 
to white dwarf and the mass loss from the system. A SSS, 
RX J0513.9-6951 (open star), whose WD mass is 1.30 Mq 
(Hachisu & Kato 2003b), locates in the initial contour. It 
should be a good candidate for the progenitor of SN la, as 
suggested in Hachisu & Kato (2001) and Hachisu & Kato 
(2003b). Three recurrent novae (U Sco, V394 CrA and 
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Fig. 1. Parameter regions producing SNe la in the 
(logP — M2) (orbital period-donor mass) plane for the WD 
-I- MS systems. The initial WD mass is I.IOMq. The 
WD -|- MS system insides the region encircled by sohd hne 
(labeled "initial") will increase its white dwarf mass up to 
1.378 Mq, where we assume a SN la explosion. The fi- 
nal state of the WD -|- MS system in the plane is encir- 
cled by dot-dashed line (labeled "final"). Filled squares indi- 
cate SN la explosions during an optically thick wind phase 
(IM2I > Mcr). Filled circles denote SN la explosions af- 
ter the wind phase, where hydrogen-shell burning is stable 
(Mcr > IM2I > |Mcr). Filled triangles denote SN la explo- 
sions after the wind phase where hydrogen-shell burning is 
mildly unstable {■^Mcr > |Af2| > gMcr). A supersoft X-ray 
sources, RX J0513. 9-6951 (open star) is plotted, whose or- 
bital period is 0.763 days (PakuU et al. 1993). Three recur- 
rent novae are indicated by a filled star for U Sco whose pe- 
riod is 1.2306 days (Schacfer & Ringwald 1995; Hachisu et al. 
2000a; Hachisu et al. 2000b), a solar symbol for CI Aq whose 
orbital period is 0.6184 days (Mennickent & Honcycutt 1995) 
and an earth symbol for V394 CrA whose orbital period is 
0.7577 days (Schaefer 1990). 
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Fig. 2. Similar to Fig. 1, but for different initial WD mass. 

CI Aql) are outside the initial contour for SNe la, while 
still inside the final region. The WD mass of U Sco was 
estimated to be 1.37 Mq (Hachisu et al. 2000a; Hachisu 
et al. 2000b), which is very close to 1.378 Mq. The MS 
mass of U Sco is 1.5 A^© (Hachisu et al. 2000a; Hachisu 
et al. 2000b), which still has enough material to trans- 
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Fig. 3. Similar to Fig. 1, but for different metallicities. 

fer onto the WD and increases the WD mass to 1.378 
Mq. So. U Sco is a very likely candidate for the pro- 
genitor of SN la (sec also Hachisu ct al. 2000a; Hachisu 
et al. 2000b). The WD masses of V394 CrA and CI 
Aql were estimated to be 1.37 Mq and 1.2 Mq, respec- 
tively. Their companion masses are still unclear. The 
best-fit companion masses for the two recurrent novae are 
1.50 Mq, while I.OMq ~ 2.OM0 is still accepted (Hachisu 
& Kato 2000; Hachisu & Kato 2003a; Hachisu, Kato & 
Schaefer 2003). So, based on our binary evolution calcu- 
lation, V394 CrA is a very likely progenitor of SN la and 
CI Aql is a possible progenitor of SN la. A further obser- 
vation to confirm their companion masses is necessary for 
judging their fates. Another very famous recurrent nova is 
V 1974 Cygni (Nova Cygni 1992), whose companion mass 
is 0.21 Mq (Deyoung & Schmidt 1994; Retter, Leibowitz 
& Ofek 1997) and orbital period is 0.0813 days (Paresce 
et al. 1995; Retter, Leibowitz & Ofek 1997). Its position 
is outside the final contour in this paper. Considering its 
low WD mass (I.OMq - I.IMq, Hachisu & Kato 2005), it 
is very unlikely for V 1974 Cygni to become a SN la (See 
Hachisu et al. (2008) for similar discussions about these 
recurrent novae). 

In figure 2, we present the initial contours for SNc la 
and the final state in (log P — M2 ) plane at the moment 
of SNe la explosion for different initial WD masses with 
Z = 0.02. This figure further confirms the likelihood for 
V394 CrA and CI Aq to be the progenitors of SNe la. 

In figure 3, we show the trend of the initial and final 
contour with metallicity. For clarity, we only show the 
case of M^Y) = I.IOMq, since the other cases give similar 
results (see Fig. 4 in Meng, Chen & Han 2009). We see 
from the figure that there is a clear trend for the initial 
and final contour to move to higher mass with metallicity. 
This is due to the correlation between stellar evolution and 
metallicity. Generally, stars with high metallicity evolve 
in a way similar to those with low metallicity but less mass 
(Umcda et al. 1999; Chen & Tout 2007; Meng, Chen & 
Han 2008). Thus, for binaries of CO WDs with partic- 
ular orbital periods, the companion mass increases with 
metallicity (see Meng, Chen & Han 2009 in details). 
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Fig. 4. Similar to figure f, but only for Case Wind. The 
initial WD mass is f.fO Mq. The big field squares show 
the initial WD -|- MS systems, while the small ones show 
their final positions at the moment of SNe la explosion. 
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Fig. 5. Similar to figure 4, but for different initial WD 
mass. Different colors represent different initial WD masses. 
The big field squares show the initial WD -|- MS systems, 
while the small ones show their final positions at the mo- 
ment of SNe la explosion. In the figure, we do not plot 
the cases with Mi < O.QMq for the clarity of appearance 
of the figure since no Case Wind occurs for the cases. 

We present the initial and final state of binary evolution 
of SNe la for Case Wind in figures 4 and 5. Figure 4 is only 
for the case of M^^} = I-IMq and figure 5 is for different 
initial WD masses. Here, we only show the case oi Z = 0.02 
as a typifier, since the other cases give similar results (see 
figure 3). In these two figure, we see that the range of 
initial MS masses with Z = 0.02 for the Case Wind is be- 
tween 2.2 Mq and 3.5 Mq. However, the range strongly 
depends on metallicity, that is: the progenitor systems 
have more massive companions for a higher Z (see figure 
3 in this paper and Fig. 4 in Meng, Chen & Han 2009). 
Hachisu et al. (2008) found that when M^j^^ < 0.9Mq, no 
Case Wind exists for Z = 0.02. We get a similar results 



No. ] 

for Z = 0.02. However, this result also rests on metallic- 
ity. For example, the mass limit for no wind is 0.8 A/© for 
Z = 0.06, while it is larger than 1.0 Mq for Z = 0.0001. 
In addition, Hachisu et al. (2008) noticed that only when 
M\ > S.OMq, WDs explode at the optically thick wind 
phase as SNe la. Their result is much different from that 
in this paper. The difference is derived from the differ- 
ent treatment on mass transfer between WD and MS and 
from different mass-loss mechanism from systems. Firstly, 
they used a analytic method for treating binary interac- 
tion. This method can not describe certain mass-transfer 
phase, in particular those occurring on a thermal time- 
scale and may overestimate mass-transfer phase on the 
thermal time-scale (Langer et al. 2000; see also the Figs. 
1 and 4 in Han & Podsiadlowski 2004). Secondly, they as- 
sume a mass-stripping effect, that is: optically thick winds 
from WD collide with secondary and strip off its surface 
layer. This effect may attenuate the mass-transfer rate be- 
tween WD and its companion, but the total mass-loss rate 
from secondary increases (Hachisu & Kato 2003a; Hachisu 
& Kato 2003b; Hachisu et al. 2008). For the two reasons 
above, more materials thus are lost from binary systems 
as the optically thick wind and stripped-off material be- 
fore supernova explosion, and then the systems are more 
likely to explode after the optically thick wind. 

The material lost as the optically think wind forms cir- 
cumstellar material (CSM) which may result in a color 
excess of SNe la (Meng et al. 2009). Since the CSM is 
very close to SN la for Case Wind, the ejecta of explosion 
may interact with the CSM, and the SN la might then be 
observed as supernovae like SN 2006X which show a vari- 
able Na I D line (Patat et al. 2007; Blondin et al. 2009). 
We will discuss this in the futiu'c paper. 

In the final contours, there is no special region for Case 
Wind (see figure 5), that is: WDs can explode as SNe 
la at any position in the final permitted region in the 
(logP — M2) plane, which is much different from the Case 
Calm and the Case Nova (see the following subsection). 

3.3. Case Calm 

If the mass-transfer rate from secondary onto WD is 
below the critical accretion rate (Equation 2) while above 
the lowest rate of steady hydrogen burning before super- 
nova explosion, i.e. M^r > IM2I > \Mcr, the WD under- 
goes steady H burning at the moment of SN la explosion. 
This is the reason why we call this "Case Calm" (Same to 
that in Hachisu et al. 2008). During the steady hydrogen 
burning phase, the WD may be observed as a SSS. The 
material lost as the optically thick wind forms CSM, but it 
has been dispersed too washy to be detected immediately 
after the SN la explosion. In figures 6 and 7, we show 
the initial and final state of binary evolution of SNe la for 
Case Calm in the (logP — M2) plane. Figure 6 is only for 
the case of = 1.1 Mq and figure 7 is for different ini- 

tial WD masses. Here, we only show the case oi Z = 0.02 
as a typifier, since the other cases give similar results (see 
also figure 3). In figure 6, it is clear that there is gap 
for both initial system and final state in the (logP — M2) 
plane, and the gap divides the Case Calm into two groups. 
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Fig. 6. Similar to figure 4, but only for Case 
Calm. The big field circles show the initial WD 
-I- MS systems, while the small ones show their fi- 
nal positions at the moment of SNe la explosion. 
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Fig. 7. Similar to Fig. 6, but for different initial 
WD masses. Different colors represent different ini- 
tial WD masses. The big field circles show the ini- 
tial WD -|- MS systems, while the small ones show 
their final positions at the moment of SNe la explosion. 

The gap in the initial contour is the region for Case wind 
and the area of the gap decrease with decreasing. 
When M^j-, < O.9M0, the gap disappear (see figure 7). 

The gap is also shown in the initial (logP — M2) plane 
of Han & Podsiadlowski (2004) and Meng, Chen & Han 
(2009), but they do not discussed it. In this paper, we 
explain its origin as foUowings. For the high-mass group 
in the initial contour of figure 6, their mass ratio is large. 
After the onset of Roche lobe overflow (RLOV), mass- 
transfer rate is so high that the system almost undergoes 
dynamically unstable. At this stage, a large amount of 
hydrogen-rich materials lose as the optically thick wind 
and hence the secondary mass decreases sharply. The 
mass-transfer rate drops after the mass ratio has been 
reversed and then the optically thick wind stops. When 
CO WD mass increases to 1.378 Mq, the secondary mass 
decreases to about 0.9 A/q, which corresponds to the low- 
mass group in the final state contour of figure 6. 
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For the low-mass group in the initial contour of figure 6, 
their mass ratio is not very large and only a small amount 
of hydrogen-rich materials lose as the optically thick wind, 
which means that most of transferred materials are accu- 
mulated on the WD. Then, the secondary mass decreases 
slightly, which forms the high-mass group in the final state 
contour of figure 6. 

In figure 7, there is a blank region in the final state re- 
gion for Case Calm (around Ml^ = 1.4A/Q-logP = 0.6), 
which means that SSSs may not be observed at this re- 
gion before SNe la explosion. This phenomena is much 
different from that of Case Wind. The blank region is 
mainly derived from the low-mass group in the initial con- 
tour. For a system belonging to the group, mass transfer 
is almost conservative, and then the period is always de- 
creasing since mass ratio is not reversed, which leads to 
the region absenting long-period systems. For a system 
in the high-mass, although its period decrease during op- 
tically thick wind, the period increases after mass-ratio 
reversion since no material loses from system and mass 
transfer becomes conservative (see the panel (h) of Fig.l 
in Han & Podsiadlowski 2004). Thus, its final period is a 
bit complicated. 

34- Case Nova 
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Fig. 8. Similar to figure 4, but only for Case Nova. 
The big field triangles show the initial WD -|- 
MS systems, while the small ones show their fi- 
nal positions at the moment of SNe la explosion. 



If mass-transfer rate from secondary onto WD is be- 
low the lowest rate of steady hydrogen burning, while still 
higher than a special value, i.e. ■^Mcr > IM2I > |-^cr, 
hydrogen shell burning is unstable to flash, and this can 
recur many times in a short period as a recurrent nova. 
We call this "Case Nova" (similar to Hachisu et al. 2008). 
During the last recurrent nova phase, the amount of the 
hydrogen-rich material lost from system is very small. On 
the other hand, the material lost as the optically thick 
wind is already much far from the center of SN la explo- 
sion. It will take 100-1000 yr for the SN la ejecta to catch 
the material lost as the optically thick wind (Hachisu et 
al. 2008), and then no signal about CSM can be detected 
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Fig. 9. Similar to figure 8, but for different initial 
WD masses. Different colors represent different ini- 
tial WD masses. The big field circles show the ini- 
tial WD + MS systems, while the small ones show 
their final positions at the moment of SNe la explosion. 



immediately after supernova explosion. In figures 8 and 
9, we show the initial and final state of binary evolution 
of SNe la for Case Nova in the (logP — M2) plane. Figure 

8 is only for the case of M^^) = l-l-^o and figure 9 is 
for different initial WD masses. Here, we only show the 
case oi Z = 0.02 as a typifier, since the other cases give 
similar results (see also figure 3). The initial region of 
Case Nova is also divided into two group, which is similar 
to that of Case Calm. The systems located in the gap 
are those for Case Wind or Case Calm. The difference 
between Case Calm and Case Nova is that the gap dis- 
appear when = O.8M0, not 0.9Mo (see figure 9). 
For a similar reason to that of Case Calm, the high-mass 
group in the initial contour forms the low-mass group in 
the final contour, while the low-mass group in the initial 
contour forms the high-mass group in the final contour. 

Similar to figure 7, there is also a blank region in figure 

9 (around Aff^ = 1.3AfQ-logP = 0.5). The origin of the 
blank region is similar to that of Case Calm. 

The recurrent nova U Sco, which is an excellent candi- 
date of SN la progenitor (Hachisu et al. 2000a; Hachisu et 
al. 2000b), is located in the middle of the final region, and 
its position in the (logP — M2) plane is just the permitted 
region for recurrent nova (see figure 9). Our model can 
explain the position of U Sco excellently. 

3.5. Companion state after SNe la explosion. 

A good way of discriminating between the many SN la 
progenitor scenarios is to search for the companion of a 
SN la in its remnant. Unless the companion is another 
WD (DD channel; in which it has been destroyed by the 
mass-transfer process itself before explosion), it survives 
and shows some special properties in its spectra, which 
is originated from the contamination of supernova ejecta 
(Marietta et al. 2000; Ruiz-Lapuente et al. 2004; Branch 
2004). Tycho's supernova, which is one of only two SNe 
la observed in our Galaxy, provides an opportunity to 
address observationally the identification of the surviving 
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Fig. 10. The final state of companion star at the mo- 
ment of SNc la explosion for different initial WD mass 
with Z = 0.02. Cross represents Tycho G, which 
is a potential candidate for the companion of Tycho's 
supernova (Ruiz-Lapuente et al. 2004; Branch 2004). 
The length of the cross represents observational error. 
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companion. Ruiz-Lapuente et al. (2004) searched the re- 
gion of the remnant of Tycho's supernova and suggested 
that Tycho G, a sun-hke star, is the companion of Tycho's 
supernova. Although, Ihara et al. (2007) argued that the 
spectrum of Tycho G do not show any special properties, 
which seems to exclude the possibility of Tycho G to be 
the companion of Tycho' supernova, the analysis of the 
chemical abundances of the Tycho G upholds the com- 
panion nature of the Tycho G (Hernandez et al. 2009). 

However, the knowledge about the companions of SNe 
la after explosions is still unclear. Generally, the super- 
nova ejecta in the single degenerate model collides into the 
envelope of its companion and strips some hydrogen-rich 
material from the surface of the companion. After the col- 
lision, the companion gains a kick velocity, which is much 
smaller than orbital velocity, and leaves explosion center 
at a velocity similar to its orbital velocity (Marietta et al. 
2000; Meng et al. 2007). 

For a companion star in the remnant of a SN la, its 



mass and space velocity can be detected directly. Since 
the stripped-off hydrogen-rich material from companion 
surface and the change of the space velocity of companion 
resulting from the collision of explosion ejecta are both 
small (Marietta et al. 2000; Mattila et al. 2005; Leonard 
2007; Meng et al. 2007), we can approximately use the 
state of a companion at the moment of supernova ex- 
plosion to present its final state after supernova explo- 
sion. In figure 10, the final states of companion stars are 
presented in Voth ~ Mf (orbital velocity-final companion 
mass) plane for different initial WD mass with Z = 0.02. 
From the figure, we can see that the range of space veloc- 
ity is from 70 km/s to 210 km/s and the mass is between 
0.6 Mq and 2.5 Mq. Tycho G just locates in this range. 
Our work can explain the observation. In the future pa- 
per, we will give detailed results about Tycho G by the 
approach of binary population synthesis. 

The effect of metallicity on the final state of the com- 
panions are shown in figure 11. In the figure, we only 
show the case of M^j-, = I.IOMq, since the other cases 
give similar results, we can see from the figure that all 
the metallicity may explain the properties of Tycho G. 
But for a detailed binary population synthesis study, only 
the case of Z = 0.02 can well explain the position of Tycho 
G in the T4rb — Mf ^ plane. Please notice our following 
paper. 

4. Discussion 

Initial and Final Contour 

In this paper, we show the initial and final parameter 
spaces for SNe la in (logP — M2) plane. The final masses 
of the companions are between 0.6 Mq and 2.5 Mq for 
Z = 0.02. However, the range of the companion mass in 
Hachisu et al. (2008) is from 1.2 Mq to about 3 Mq. 
This difference is directly derived from the different ini- 
tial contours, which is mainly determined by the different 
treatment of thermal timescale mass transfer and different 
mass-loss mechanism. They use a simply analytic method 
to estimate the mass-transfer rate, which leads to a differ- 
ent initial contour for SNe la, especially for low-mass WD 
(Han & Podsiadlowski 2004). In addition, they assumed 
a mass-stripping effect, and then, the initial companion 
mass may be as large as 8 A/q, which depends on the 
efficiency of the mass-stripping effect. Since their final 
mass of companion is always larger than 1.2 A/q, it is dif- 
ficult for their model to explain the properties of Tycho 
G; which remains consistent with the surviving compan- 
ion of Tycho's supernova (1 Mq, Ruiz-Lapuente et al. 
2004; Hernandez et al. 2009). Since our results are based 
on detailed binary evolution calculation with latest input 
physics, our results are more realistic and the final con- 
tours in this paper are more likely to approach the real 
ones than those in Hachisu et al. (2008) (see also in Han 
& Podsiadlowski 2004) . Because of the difference range of 
final companion as mentioned above, Tycho G can be a 
likely progenitor of Tycho's supernova in our model, but 
can not be in the model of Hachisu et al. (2008). 
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4-2. Supernova like SN 2002ic and V Sge 
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Fig. 12. Initial region for SNc la with different metallici- 
ties. Here, only the cases of Af^^ = 1.2OM0 are shown. 
The black solid line shows the initial region of the progen- 
itor for supernovae like SN 2002ic. Two supersoft X-ray 
sources, RX J0513. 9-6951 (open star) and V Sge (filled star), 
are plotted, whose masses are 2.7 Mq and 3.5 Mq (Hachisu 
& Kato 2003b; Hachisu &; Kato 2003c), and whose orbital 
period are 0.76 days and 0.51 days, respectively (PakuU 
et al. 1993; Herbig et al. 1965; Patterson ct al. 1998). 

Until SN 2002ic was discovered (Hamuy et al. 2003), 
it was long believed that there are no hydrogen lines in 
the spectra of SNe la. The strong hydrogen lines in the 
spectra of SN 2002ic were explained by the interaction be- 
tween the SN ejecta and the circumstellar material (CSM) 
(Hamuy et al. 2003). Recently, two co-twins of SN 2002ic 
(SN 2005gj and SN 2006gy) were also found (Aldering 
et al. 2006; Ofck et al. 2007). To explain these rare 
objects, many models were suggested and here we just 
list some of them: Hachisu et al. (1999b); Hamuy et al. 
(2003); Livio & Riess (2003); Chugai & Yungelson (2004); 
Han & Podsiadlowski (2006); Hachisu ct al. (2008); Lii 
et al. (2009). Among all the models listed above, the 
delayed dynamical instability model suggested by Han & 
Podsiadlowski (2006) is more interesting and some pre- 
dictions from the model seem to be consistent with obser- 
vations, especially in a sense of the birth rate and delay 
time of these rare objects (Aldering ct al. 2006; Prieto 
et al. 2007). In the scenario of Han & Podsiadlowski 
(2006), SN 2002ic might be from the WD -f- MS channel, 
where the CO WD accretes material from its relatively 
massive companion ('^ 3.OM0), and increases its mass to 
^ 1.3OA'/0 before experiencing a delayed dynamical insta- 
bility. Adopting their model, Meng, Chen & Han (2009) 
assumed that the CO WD can increase its mass to 1.378 
Mq and explode as a SN 2002ic-like case if the mass of 
the CO WD exceeds 1.30 Mq before the delayed dynami- 
cal instability, and then predicted that supernovae like SN 
2002ic may not be found in extremely low-metallicity envi- 
ronments. Under the same assumption to Meng, Chen & 
Han (2009) , we show the permitted region for the delayed 
dynamical instability model with various metallicities in 
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figure 12. For the simplicity of the assumption, the real 
region for the delayed dynamical instability model may 
be larger than that shown in figure 12 (see Fig. 5 in 
Han & Podsiadlowski 2006 and Fig. 2. in Meng, Chen 
& Han 2009). It is clear that the permitted region for 
supernova like SN 2002ic locates above the contours for 
normal supernova. There is a overlap between supernova 
like SN 2002ic and normal SN la, since a progenitor of 
supernova like SN 2002ic with low metallicity may have 
a lower mass than that of normal supernova. Whatever, 
more evidence is necessary to confirm the scenario in Han 
& Podsiadlowski (2006), especially to find a progenitor 
system at present. 

V Sge is a well observed quasi-periodic transient SSS 
in our Galaxy. It switches optical state between a high 
iy ~ 11 mag) and low {y ^ 12 mag) state, and during 
the low state, very soft and very weak X-ray can be de- 
tected (Greiner & van Tecschng 1998). The mass- loss 
rate for V Sge can be as large as ~ lO~^M0yr~^ indi- 
cated by radio observation. (Lockley et al. 1997; Lockley 
et al. 1999). Hachisu & Kato (2003c) use optically thick 
wind to explain the observational property of V Sge excel- 
lently, and in their model, the mass-loss rate may reaches 
as high as ^ lO~^M0yr~^, which is consistent with radio 
observations. Hachisu & Kato (2003c) also estimated the 
WD mass as Mwd ~ 1.25M0 and its companion mass as 
M2 ^ 3.5M0, and then they suggested that V Sge will ex- 
plode as a SN la after about 10~^yr. However, they must 
assume a mass-stripping effect, otherwise the position of 
V Sge in (logP — M2) plane will beyond the permitted re- 
gion for stable mass transfer (Hachisu et al. 2008). We get 
a similar result to that in Hachisu et al. (2008) that V Sge 
locates in a forbidden zone for stable mass transfer if there 
is no special assumption (see in figure 12). However, it lo- 
cates at the boundary of the contour for supernovae like 
SN 2002ic as indicated in Han & Podsiadlowski (2006). 
Considering the smaller region here than a real one and 
a moderately larger WD mass of V Sge than that in this 
paper, V Sge is a potential progenitor system that will 
undergo the delayed dynamical instability as suggested 
in Han & Podsiadlowski (2006). V Sge thus would not 
explode as a normal SN la, but one like SN 2002ic, if 
the model in Han & Podsiadlowski (2006) is appropriate. 
Then, we might regard "V Sge-type star" as the progen- 
itor of supernovae like SN 2002ic if it belongs to WD + 
MS system. 

The lifetime of V Sge-type star is typically ^ lO^yr, 
which is mainly due to the time-averaged wind mass- 
loss rate of ~ IQ-^Mqyt-'^ (Han & Podsiadlowski 2006). 
Since no more than 1 in 100 SNe la belong to the sub- 
group of 2002ic-like supernova (Han & Podsiadlowski 
2006; Aldering et al. 2006; Prieto et al. 2007), there should 
be several V Sge-type star belonging to WD -I- MS sys- 
tem in our Galaxy if we take the Galactic birth rate of 
SNe la as 3-4x 10~^yr~'^ (van den Bergh & Tammann 
1991; CappcUaro & Turatto 1997). Steiner & Diaz (1998) 
listed four V Sge-type stars in the Galaxy and discussed 
their similar spectroscopic and photometric properties. At 
present, their companion masses have not been identified. 
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while their orbital periods is in the range of 0.2-0.5 days, 
which falls into the orbital periods predicted by our model 
(see the final regions in figure 3). Please notice that not 
all SSSs belong to WD + MS systems and the nature of 
the compact star still remains an open question (Steiner 
& Diaz 1998). White dwarfs, neutron stars, black holes, 
or even carbon main-sequence stars are all possible candi- 
dates (Steiner & Diaz 1998). Another alternative for the 
SSSs is WD -I- He star system, which may also contribute 
to a few part of SNe la (Wang et al. 2009). Thus, consid- 
ering the uncertainty of the nature of the four V Sge-type 
star, our estimation on the number of V Sge-type stars is 
consistent with observations. 

However, although the life, the mass-lose rate and the 
number of V Sge-type star seem to match with the pre- 
diction from the delayed dynamical instability model, the 
case of = 1.2Mq is the only one in our model grids 

to account for the position of V Sge in the (logP — M2) 
plane, while there exist arguments on whether the mass of 
a CO WD may be as large as 1.2Mo and whether the CO 
WD of I.2M0 may explode as a SN la (Umeda et al. 1999; 
Meng, Chen & Han 2008). So, obtaining a conclusion on 
whether there is a relation between V Sge-type star and 
the delayed dynamical instability model is premature, and 
V Sge is only a potential candidate for the progenitor of 
supernovae like SN 2002ic. 

5. Summary and Conclusions 

Adopting the prescription in Hachisu et al. (1999a) for 
the mass accretion of CO WDs, we have carried out de- 
tailed binary evolution calculations for the progenitors of 
SNe la with different metallicity in the single-degenerate 
channel (the WD + MS channel) and obtained the initial 
and final parameters in the (logP — A/2) plane. Our model 
may explain the positions of some famous recurrent novae 
in the (logP - M2) plane, as well as a SSS, RX J0513.9- 
6951. Our model can also explain the space velocity and 
mass of Tycho G, which is now suggested to be the poten- 
tial companion star of Tycho's supernova (Ruiz-Lapuente 
et al. 2004; Hernandez et al. 2009). Based on the de- 
layed dynamically instable model in Han & Podsiadlowski 
(2006), we might regard V Sge as a potential progenitor 
of supernova like SN 2002ic. 
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